Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN M?eX style file vl.4) 



The history of star formation in dusty galaxies 



A. W. Blain/ Ian Smail,^ R. J. Ivison^'^ and J. -P. Kneib^ 

^ Cavendish Laboratory, Madingley Road, Cambridge, CB3 OHE, UK. 

^ Department of Physics, University of Durham, South Road, Durham, DHl 3LE, UK. 

^ Institute for Astronomy, Department of Physics & Astronomy, University of Edinburgh, Blackford Hill, Edinburgh, EH9 3HJ, UK. 
^ Department of Physics & Astronomy, University College London, Cower Street, London, WCIE 6BT, UK. 
^ Observatoire Midi- Pyrenees, 14 Avenue E. Belin, 31400 Toulouse, France. 



1 February 2008 



ABSTRACT 

A population of distant dusty galaxies emitting in the submillimetre waveband 
has recently been detected using the Submillimetre Common-User Bolometer Array 
(SCUBA) camera on the James Clerk Maxwell Telescope (JCMT). This population 
can be used to trace the amount of high-redshift star-formation activity that is ob- 
scured from view in the optical waveband by dust, and so is missing from existing 
inventories of star formation in the distant Universe. By including this population 
we can construct a complete and consistent picture of the history of star formation. 
The evolution of obscured star formation at redshifts less than unity is constrained 
by mid- and far-infrared counts of dusty galaxies. Activity increases with redshift z 
as (1 -|- z)'' with 7 ~ 4, consistent with the form of evolution found in the optical 
waveband by the Canada-France Redshift Survey (CFRS) to z < 1. The form of evo- 
lution at higher redshifts is constrained by both faint SCUBA counts and the intensity 
of background radiation in the millimetre/submillimetre waveband. We find that the 
total amount of energy emitted by dusty galaxies is about four times greater than 
that inferred from restframe ultraviolet observations, and that a larger fraction of this 
energy is emitted at high redshifts. The simplest explanation for these results is that 
a large population of luminous strongly-obscured sources at redshifts z < 5 is missing 
from optical surveys. We discuss the possible contribution of obscured active galactic 
nuclei to the submillimetre-wave background and counts. More accurate constraints 
on the history of star formation will be provided by determinations of the counts in 
several submillimetre wavebands and crucially by a reliable redshift distribution of the 
detected galaxies. 

Key words: galaxies: evolution - galaxies: formation - cosmology: observations - 
cosmology: theory - diffuse radiation - radio continuum: galaxies 



1 INTRODUCTION 

A history of star formation in galaxies over roughly the last 
half of the Hubble time has been inferred with reasonable 
accuracy by combining the star-formation rate (SFR) in the 
local universe (Gallego et al. 1996; Gronwall 1998; Tresse 
& Maddox 1997; Treyer et al. 1998), the results of spectro- 
scopic redshift surveys out to z ~ 1 (Lilly et al. 1996) and the 
results of deep photometric surveys, such as the Hubble Deep 
Field, in which Lyman-dropout selection techniques (Steidel 
et al. 1996a) are used to determine photometric redshifts: 
see Madau et al. (1996) and Madau, Delia Valle & Panagia 
(1998). This history of unobscured star formation in galaxies 
can be described adequately by a power-law increase in the 
comoving density of star-formation activity with redshift z 
of the form (1 + z)'' with 7^4 out to z ~ 1, followed by a 



decrease at higher redshifts. Observations made in the near- 
infrared waveband by Connolly et al. (1997) have provided 
further support for the suggestion that the SFR peaks at 
z ~ 1-2. This standard picture of the star-formation history 
is shown in Fig. 1. It accords with the form of evolution of 
the luminosity density inferred from observations of radio 
sources (Cram 1998; Dunlop 1998) and quasars (Boyle & 
Terlevich 1998). 

This work has greatly increased our knowledge of the 
evolution of star formation activity in the Universe; how- 
ever, the points plotted in Fig. 1 in fact only represent firm 
lower limits to the true SFR at each epoch. Extinction due 
to any interstellar dust in the optically selected observed 
galaxies suppresses their flux densities in the rest-frame op- 
tical and ultraviolet wavebands, and so leads to an under- 
estimation of their SFR, sometimes by quite considerable 
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factors 10). Using estimates of obscuration in Lyrnan- 
dropout galajcios derived from restframe ultraviolet spectra, 
Pettini et al. (1998a,b) have recently revised the global SFR 
at redshifts greater than 2 upwards by a factor of about 5, 
as shown by the open boxes in Fig. 1. However, it is difflcult 
to correct optical and ultraviolet spectra reliably for the ef- 
fects of dust, and so these revisions are rather uncertain. 
In contrast, the associated rest-frame mid- and far-infrared 
thermal emission from the interstellar dust grains that are 
heated by absorbing optical and ultraviolet photons can bo 
detected directly and used to place limits on the dust mass 
and SFR of the galaxies. 

The first observations of dust emission from external 
galaxies were made by the IRAS satellite. Nearby galax- 
ies selected on the grounds of blue colours in the optical 
waveband were found to emit typically 60 per cent of their 
total luminosity in the far-infrared waveband (Mazzarella & 
Balzano 1986). The redshift distribution of the population of 
galaxies detected at a wavelength of 60 /im by IRAS extends 
out to z ~ 0.2, and is consistent with pure luminosity evo- 
lution of the form (1 -|- z)^'^ (Saunders et al. 1990; Oliver, 
Rowan- Robinson & Saunders 1992). More recent observa- 
tions at wavelengths of 15 ^m and 175 ^m using the more 
sensitive ISO satellite (Kawara et al. 1997; Rowan-Robinson 
et al. 1997; Hammer & Flores 1998; Lagache et al. 1998) have 
confirmed that strong evolution continues out to « ~ 1, a re- 
sult that is broadly consistent with the form of evolution de- 
rived from optical observations (Lilly et al. 1996). However, 
although they trace the star-formation history over at least 
half the Hubble time, these observations do not strongly 
constrain the SFR in the early Universe, an era crucial for 
distinguishing between competing galaxy formation models 
(Cole et al. 1994). One of the most promising routes to probe 
this regime are deep observations at still longer wavelengths, 
in the millimetre and submillimetre wavebands. 

Early bounds on the intensity of background radiation 
in the submillimetre and far-infrared wavebands from the 
FIRAS instrument on the COBE satellite (Wright et al. 
1994; Fbcsen et al. 1996) provided limits on the amount and 
epoch of obscured star formation activity in an hierarchi- 
cal model (Blain & Longair 1993b) . Much more information 
about the role of dust was provided by Puget et al. (1996), 
who detected a considerably larger isotropic non-Galactic 
background radiation signal in the same data by using a dif- 
ferent model of Galactic emission. The results indicate that 
a large amount of far-infrared emission must be produced 
by dust at z ^ 1 (Burigana et al. 1997). Recently, this early 
estimate of the background radiation intensity has been in 
large part confirmed by independent analyses (Guiderdoni 
et al. 1997; Fixson et al. 1998), and by complimentary data 
obtained using the DIRBE instrument (Hauser et al. 1998; 
Schlegel, Finkbeiner & Davis 1998). A significant amount 
of star formation activity at redshifts greater than unity, at 
least in dense environments, is apparently required to ex- 
plain observations of metal enrichment in intracluster gas 
at z ~ 0.5 (Mushotzky & Loewenstein 1997). 

The population of submillimetre-luminous galaxies that 
produce the far-infrared/submillimetre-wave background 
signal were detected by Smail, Ivison & Blain (1997, here- 
after SIB). The counts of faint galaxies in the submillimetre 
waveband are very sensitive to the form of galaxy evolu- 
tion at high redshift because the fe-corrections expected for 
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Figure 1. The star-formation rate (SFR) as a function of redshift 
tsfrX-z), as inferred from ultraviolet /optieal/near-infrared obser- 
vations by Gallego et al. (1996; filled triangle), Gronwall (1998; 
diagonal cross), Treyer et al. (1998; open triangle), Tresse & Mad- 
dox (1998; empty circle), Lilly et al. (1996; stars), Connolly et 
al. (1997; filled squares) and Madau et al. (1996; filled circles). A 
modified high-redshift SFR, inferred after a correction for dust ex- 
tinction by Pettini et al. (1998a), is shown by the empty squares. 
The thick dotted curve is the SFR inferred from radio data by 
Dunlop (1998). The evolution of the luminosity density of AGN 
at a wavelength of 280 nm (Boyle & Terlevich 1998) is shown by 
the thin dotted line; it is converted into an equivalent SFR, rgpf^ 
assuming that tsfr = 10"'^ Mq yr~^ Mpc^'^ is equivalent to a 
280-nm luminosity density of 1.7 x lO" WHz"! Mpc-^. 

distant galaxies are very large and negative. Hence high- 
redshift galaxies are likely to make a substantial contribu- 
tion to the counts (Blain & Longair 1993a), which in turn 
provide a sensitive test of galaxy evolution models in the 
early Universe. 

SIB exploited the especially strong magnification bias 
expected in rich clusters of galaxies in the submillimetre 
waveband (Blain 1997) to increase the detoctability of the 
population of background galaxies. Unless the cluster galax- 
ies have extremely high and undetected SFRs, the contribu- 
tion of flux densities from galaxies within the lensing clusters 
should be at least an order of magnitude less than that from 
background sources. The observations were carried out using 
the sensitive SCUBA camera at the James Clerk Maxwell 
Telescope (Holland et al. 1999). The results of the completed 
cluster survey and its follow up are described elsewhere (Ivi- 
son et al. 1998; Smail et al. 1998; Blain et al. submitted; 
Smail et al. in preparation). The first results of other blank- 
field surveys (Barger et al. 1998; Bales et al. 1998; Holland 
et al. 1998b; Hughes et al. 1998) are now available. 

We use those 850-^m observations and other sub- 
millimetrc/far-infrarcd observations at wavelengths of 
175 //m (Kawara ct al. 1997) and 2.8 mm (Wilner & 
Wright 1997), the intensity of background radiation in the 
millimetre/submillimetre/far-infrared wavebands, and ex- 
isting information about the evolution of dusty galaxies at 
low redshifts derived from 60-yum observations (Oliver et al. 
1992) to investigate the history of star-formation activity in 
the Universe that is obscured from view in the optical wave- 
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band by dust. Throughout the paper we use the comoving 
density of star-formation activity to represent the luminosity 
density of the Universe, in order that the results of obser- 
vations in the optical and submillimetre wavebands can be 
readily compared. However, non-stellar omission from active 
galactic nuclei (AGN) contribute to the total emitted radia- 
tion in the submillimetre waveband. This issue is discussed 
in more detail in Sections. 

We review the relevant observations and existing lim- 
its in Section 2, outline our strategy for determining the 
star formation history in Sections, and discuss the self- 
consistency of the results in Section 4. In Sections we 
present the resulting star formation histories and review 
them in the context of optical observations. In Section 6 
we investigate the most promising routes to obtaining bet- 
ter constraints on the history of star-formation activity in 
the millimetre/submillimetre and far-infrared wavebands: 
sec also Blain (1998c). A value of Hubblc's constant Ho = 
50kms~^ Mpc~^, a density parameter Oq = 1 and a cosmo- 
logical constant = are assumed. 



2 OBSERVATIONS 

2.1 Mid/far-infrared source counts 

The whole sky was observed at a wavelength of 60 /im in 
the IRAS survey, and counts of galaxies were determined to 
various flux density limits (Hacking & Houck 1987; Rowan- 
Robinson et al. 1990; Saunders et al. 1990; Gregorich et al. 
1995; Berlin, Dennefcld & Moshir 1997). At this wavelength 
the emission from galaxies is dominated by the thermal ra- 
diation of dust grains, heated by absorbing the optical and 
ultraviolet light from both young high-mass stars and AGN. 
The most luminous class of sources detected by IRAS appear 
to be powered by stars and AGN in the ratio 2:1 (Sanders 
& Mirabel 1996). The redsliift distribution of the faintest 
galaxies in the survey extends out to a redshift of about 
0.2. Saunders et al. (1990) and Oliver et al. (1992) found 
that pure luminosity evolution of the form (l-|-z)^'^ or pure 
density evolution of the form (1 -|- z)^'^ can account for the 
results. The inferred form of luminosity evolution is compa- 
rable with that determined for radio-galaxies, quasars and 
optically-selected star-forming galaxies (see e.g. Dunlop & 
Peacock 1990; Hewett, Foltz & Chaffee 1993; Lilly et al. 
1996; Boyle & Terlcvich 1998 and Dunlop 1998). However, 
it is possible to fit the counts adequately by a range of evo- 
lution parameters if the typical dust temperature is varied: 
see Section 3.3. 

The mid- and far-infrared counts of more distant galax- 
ies are being determined from the results of smvcys using 
the ISO satellite; for example the observations of Kawara 
et al. (1997) and Lagache et al. (1998) at a wavelength of 
175 /im. At shorter wavelengths of 7 and 15 /im, Oliver et 
al. (1997) and Hammer & Flores (1998) derived counts from 
ISO observations of the Hubble Deep Field and CFRS fields 
respectively. However, it is difficult to translate between such 
a measurement of hot dust and the cooler dust that dom- 
inates the energy emission of a galaxy, and so the inferred 
global SFR (Rowan-Robinson et al. 1997) is very uncertain 
(Hughes 1996). Spectral features in dust emission in the mid- 
infrared waveband are also likely to affect the results (Xu 



Wavelength / jum 
2000 1000 500 200 100 




300 500 1000 3000 

Frequency / GHz 



Figure 2. The intensity of background radiation in the millime- 
tre, submillimetre and far-infrared wavebands, as inferred by: 
Pugct ct al. (1996) - thin solid lines with error bars at the ends; 
Fixsen ct al. (1998) - dotted hnc: Schlcgcl et al. (1998) - stars; 
Hauser et al. (1998) and Dwek ct al. (1998) - thick solid crosses. 
The filled circles represent lower limits to the background inten- 
sity inferred from the counts of SIB. The connected arrows point 
to the corresponding background intensity that is extrapolated 
using the count model described in SIB and Blain, Ivison & Smail 
(1998a, hereafter BIS). 

et al. 1998). Much more information about the population 
of mid-infrared-luminous galaxies will be proved in 1999 by 
the Wide-Field Infrared Explorer (WIRE) satellite (Shupe 
et al. 1998). 

2.2 The intensity of background radiation 

The intensity of diffuse background radiation in the mil- 
limetre, submillimetre and far-infrared wavebands has been 
constrained by observations of the whole sky using the FI- 
RAS (Fixsen ct al. 1994) and DIRBE (Hauser ct al. 1996) 
instruments aboard the COBE satellite (Puget et al. 1996; 
Guiderdoni et al. 1997; Fixsen et al. 1998; Hauser et al. 
1998; Schlegel et al. 1998). At wavelengths longer than 
about 400 iMm the cosmic microwave background radiation 
(CMBR) dominates the spectrum detected by FIRAS. After 
subtracting the CMBR signal and the emission from Galac- 
tic dust, the isotropic signal from extragalactic sources re- 
mains. This subtraction is a complex process, and will im- 
prove in the future. The Planck Surveyor satellite (Bersanelli 
et al. 1996) will also obtain much more accurate maps of 
the submillimetre-wave sky. However, our knowledge of the 
intensity of extragalactic background radiation in the mil- 
limetre/submillimetre and far- infrared wavebands is now 
substantially complete (see Fig. 2). Consistent limits to the 
background radiation intensity at shorter wavelengths have 
been determined by Dwek & Slavin (1994), Biller et al. 
(1998) and Dwek et al. (1998). The energy-normalized back- 
ground spectrum i//^ hence has a peak at a wavelength be- 
tween about 100 and 200 lira.. 
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2.3 Millimetre/submillimetre-wave source counts 

Millimetre/submillimetre-wave observations using both 
single-antenna telescopes and interferometer arrays can be 
used to detect and limit the counts of discrete sources. At 
a wavelength of 2.8 mm, Wilner & Wright (1997) used the 
BIMA interferometer to imago the Hubble Deep Field to a 
la sensitivity of 0.71 mjy. No sources were detected at a sig- 
nificance greater than 5a, which translates into a 3a upper 
limit to the counts of 160 deg~^ at a flux density of 2 mJy. 

At 850 iim SIB detected six sources, which corresponds 
to a surface density of (2.5 ± 1.4) x 10'' dog^^ at a flux den- 
sity of 4mJy. The brightest 850-/xm source (Ivison et al. 
1998b) was also detected at a wavelength of 450 /im, indi- 
cating a very wealc 3cr limit of about 10^ deg~^ to the 450- 
/Ltm counts at a flux density of 80 mJy. Observations of five 
additional clusters (Smail ct al. in preparation) have now 
extended the area of this survey to greater than 0.01 deg^. 
The surface density of discrete submillimetre-wave sources 
detected in the complete survey is fully consistent with mea- 
surements of the background radiation intensity discussed 
above (SIB). In deep images of dusty disks around stars, 
Holland et al. (1998) detected four unknown, probably ex- 
tragalactic objects, indicating a count of (10±6) x 10^ deg~^ 
at a flux density of 8 mJy. At shorter wavelengths, Phillipps 
(1997) presented a 350-/itm 8.6-arcmin^ map made using the 
SHARC bolometer array at the Caltech Submillimetre Ob- 
servatory, which contains no signiflcant sources and reaches 
a la sensitivity of about 100 mJy. 

The results of blank-fleld surveys in the Hubble Deep 
Field (Hughes ct al. 1998), CFRS fields (Bales et al. 1998) 
and Hawaii Deep Survey fields (Barger et al. 1998), and 
independent observations of clusters will reduce the errors 
on the 850-/L(m counts (see Fig. 11c for more details). The 
large complete cluster survey of Smail ct al. (in prepara- 
tion) will probably yield the most accurate constraints; the 
source counts derived from this large sample are discussed 
in a forthcoming paper (Blain et al. submitted). 

2.4 Early results 

Measurements of the far-infrared background radiation in- 
tensity, and deep 750 observations in the mid-infrared wave- 
band have led to estimates of the star- formation history that 
is obscured by dust, as shown in Fig. 3. The comoving lumi- 
nosity density emitted by dust in the far-infrared waveband 
ppiR can be related to the density of metals generated in 
high- mass stars by the present epoch f2m(0) in a popula- 
tion of galaxies that undergo their first generation of star 
formation at a redshift zq, by, 

pp:a(.) = 0.007pc^g^(l + .)^/^ (1) 

in an hierarchical model (Blain & Longair 1993b). The 
smoothed density of the Universe is p. The relationship is 
similar in other models. The background radiation intensity 
shown in Fig. 2 can be explained by the generation of about 
fim = 10~* since a redshift zo = 5 in this model. The esti- 
mated star-formation rate, assuming a standard conversion 
factor between far-infrared luminosity and the SFR (see Sec- 
tion 4) provides the star formation history shown in Fig. 3. 
The limits derived in a similar analysis by Burigana et al. 
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Figure 3. Histories of star formation derived from observations 
of dust emission and metal abundance In clusters. The data points 
from Fig. 1 arc rcplottcd for comparison. Blain & Longair (1993b) 
— dashed line - and Burigana ct al. (1997) - solid and dotted 
lines with arrows — imposed limits from background radiation 
intensities. The SFR inferred by Rowan-Robinson et al. (1997) 
and Hammer & Flores (1998) based on ISO observations at 7 
and 15 /im are shown by the solid crosses and empty stars re- 
spectively. The star-formation history associated with the count 
model discussed in SIB and BIS Is shown by the thick solid line. 
The dot-dashed line shows the estimate derived by Mushotzky & 
Loewenstein (1997) based on the metalllcity of lntrax;luster gas. 

(1997) are also shown. Counts of galax;ies at 7 and 15 /^m 
can be used to infer a star-formation rate at a ^ 1 (Rowan- 
Robinson ct al. 1997; Hammer & Flores 1998), although po- 
tentially with a very large uncertainty (see Fig. 3). 



3 DETERMINING THE HISTORY OF 
OBSCURED STAR FORMATION 

3.1 Strategy 

Even before the completion of redshift surveys of 
submillimetre-selected galaxies we can still obtain useful 
constraints on the star-formation history at high redshifts 
by combining observations of discrete sources with measure- 
ments of the intensity of background radiation in the mil- 
limetre, submillimetre and far-infrared wavebands to derive 
self-consistent models for the evolution of the population 
of dusty galaxies. At very low redshifts, z < 0.2, these are 
constrained to match the normalization and evolution of the 
60-/im counts of IRAS galax.ies. Note that redshifts less than 
two are described as 'low' in this paper. 

First, an estimate of the typical dust temperature and 
the slope of the evolution function at low redshifts is ob- 
tained by fitting a pair of functional forms, which asymptote 
to power laws in (1 -|- 2) at low redshift, to both the 60- and 
175-/xm counts. Most of the degeneracy between the values 
of dust temperature and evolution parameter is broken by 
including information from the 175-/Lim counts, and so a ro- 
bust estimate of the form of evolution of dusty galaxies at 
low redshifts can be obtained. 
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Secondly, the resulting models of the evolving popula- 
tion of dusty galaxies are extrapolated out to high redshifts, 
and are constrained to fit the background radiation spec- 
trum and submillimetre-wave counts. 

Thirdly, we employ a simple chemical evolution model 
to check whether or not the best-fitting models are self- 
consistent. A sufficiently large mass of dust must be gener- 
ated in each model to absorb early starlight and to account 
for the large amounts of far-infrared emission (Eales & Ed- 
munds 1996, 1997; Fraycr & Brown 1997), subject to the 
constraint on the temperature of emitting dust derived at 
low redshifts. The total mass of heavy elements produced 
in high-mass stars throughout the history of the Universe 
(Savage & Scmbach 1996) is also checked for consistency 
with observations. 

In order to avoid introducing a large number of free 
parameters, we aim to explain all the data using a single 
population of dusty galaxies, described by the 60-pm lumi- 
nosity function, with a fixed dust temperature Td and cmis- 
sivity fed- Saunders ot al. (1990) fit the GO-jim counts using 
both a warm and cool component, and Pearson & Rowan- 
Robinson (1996) adopt a similar scheme of cool 'cirrus' and 
warm 'starburst' galaxies. Although distinct populations of 
dusty galaxies exist in the local Universe (Sanders & Mirabel 
1996), in this paper we are investigating the global evolu- 
tion of galaxies, driven by observations that are sensitive 
to galaxies out to redshifts of about 10. Hence, while such 
multi-component models may be required to explain more 
detailed observations in the future, we do not believe that 
the introduction of additional parameters is justified by the 
data at present. Here we fit four observations - the evolu- 
tion of the 60-/jm count, the 175- and 850-/jm counts and 
the background radiation spectrum - using a model that 
contains four parameters - a typical dust temperature Td, a 
low-redshift slope of the evolution function p, and two red- 
shifts, 2:max and Z(], the redshifts beyond which evolution 
slows down and at which the first galaxies appear respec- 
tively. 

3.2 Spectral energy distribution 

We assume that the frequency-dependent emissivity func- 
tion of dust ku can be described in the millimetre and sub- 
millimetre wavebands by a simple power-law function v^'^ 
(Hughes 1996; Ivison et al. 1998a). The dust emission spec- 
trum is hence a modified Planck function with a Rayleigh- 
Jeans power-law index of 3.5. At shorter far-infrared wave- 
lengths this function would decline exponentially as a func- 
tion of frequency, turning over at a wavelength of 130±50 /im 
for a typical dust temperature of order 50 ± 20 K (Hughes 
1996). A power-law spectrum of the form is grafted onto 
the modified Planck function at wavelengths shorter than 
that for which the equivalent power-law index of the modi- 
fied Planck function is less than -a. The slope of this short- 
wavelength spectrum reflects the fraction of dust grains in 
the model galaxy at temperatures larger than the defining 
temperature in the Planck function. A reasonable fit to ra- 
diative transfer models of dusty star-forming galaxies is pro- 
vided by a value a — 2.2. These models arc important for 
interpreting observations at wavelengths shorter than 50 /im 
but not very relevant here (see Goldschmidt et al. 1997 and 
Xu et al. 1998). 



3.3 Fitting the low-redshift counts 

The 60-^m counts of /idyls' galaxies can be used to constrain 
the form of evolution at redshifts z <^ 1 (Saunders et al. 
1990; Oliver et al. 1992). If the luminosity function ^{L,z) 
is parametrized using the two functions n{z) and g{z), which 
describe density and luminosity evolution respectively, then, 



$(!/, 0) = n(2)$ 



,0 



(2) 



Both n{z) = (1 -|- z)^'"^ and g{z) = 1 and n{z) = 1 and 
g{z) = (1 -|- z)^'^ were found to provide good fits to the 
counts. The first of these cases is pure density evolution, 
and the second is pure luminosity evolution. Only models of 
pure luminosity evolution are considered here, because it is 
impossible to reproduce the counts of distant submillimetre- 
selected objects in a model of pure density evolution, with- 
out overproducing the intensity of far-infrared background 
radiation and the density of heavy elements generated in nu- 
cleosynthesis by a huge amount: between 50 and 100 times. 
This is because the effects of pure density and pure lumi- 
nosity evolution are almost the same in the calculation of 
the background intensity and heavy element abundance, but 
luminosity evolution has a more significant effect in the cal- 
culation of the counts. This is not to say that no density 
evolution occurs, but luminosity evolution must be domi- 
nant if the existing data is accurate. 

Two classes of evolution models are considered: an 
'anvil' model defined by the evolution function. 



9a.{z) = 



{I + Zf, Z<Zn 
(1 + 2max)^, Z > 



and a 'peak' model, defined by the function. 



g^{z) = 2 ( 



1 + cxp ■ 



(1 z)P+(2^°-''' \ 



(3) 



(4) 



In both models no galaxies exist beyond a redshift of zo- At 
very low redshifts gp has the asymptotic form (1 -|- zY- 

These functional forms allow rather general and plau- 
sible evolution functions to be simply described: compare 
Cole et al. (1994). The peak model can be used to de- 
scribe the star-formation history derived from optical and 
near-infrared observations in Fig. 1 adequately: see Table 1 
and Fig. 8. By requiring an evolving IRAS 60-/Lim luminos- 
ity function to fit both the form of evolution of the 60-/jm 
counts and 175- /nm ISO count, best-fitting values of the dust 
temperature Td and low-redshift evolution parameter p can 
be determined, as shown in Fig. 4. The derived best-fitting 
value of the dust temperature Td = 38 ± 4 K in both the 
anvil and peak models. In the anvil model the best-fitting 
evolution parameter p — 3.8±0.2, and in the peak model the 
best-fitting p = 3.9 ± 0.2. A slightly stronger form of evolu- 
tion is required in the peak model because the exponential 
term reduces the strength of the evolution at redshifts less 
than Zmax- Note that these slopes are very similar to the 
value required to fit the rate of evolution of optically se- 
lected star-forming galaxies determined at a wavelength of 
280 nm by Lilly et al. (1996): they found p = 3.9±0.75. The 
estimated dust temperature appears to be consistent with 
the typical value of less than 40 K derived for a sample of 
very distant quasars by Benford et al. (1998). The 60-/im 
counts predicted by the best-fitting models are compared 
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Figure 4. The results of fitting the IRAS 60-/im counts and the ISO 175-fim count (Kawara ct al. 1997) using the anvil and peak 
models: equations (3) and(4) respectively. In (a) and (b) the solid contours show values of x^/N that are 1, 2, 3, 4 and 5a away 
from the local minimum values derived from the IRAS count data; the thick and thin dashed lines trace the locus of parameters 
for which the exact, half and twice the 750 count is reproduced respectively. Likelihood contours, representing the same values, 
obtained by fitting the models jointly to the IRAS and 750 counts are shown in (c). The best-fitting parameters in the anvil and 
peak model are represented by a solid and empty circle and correspond to x^ values of 2.63 and 2.73 respectively. As the IRAS and 
ISO sources are expected to be at low redshifts, there is little difference between the predictions of the models. 



with the IRAS data and the predictions of a non-evolving 
model in Fig. 5. 

Equally good fits to low-redshift evolution of dusty 
galaxies are provided by the peak and anvil models. How- 
ever, the 60-/Ltm counts provide no significant restrictions 
on the form of evolution at 2 ^ 1, and the 750 counts are 
not expected to include a significant proportion of galaxies 
at redshifts greater than 2. Hence, the submillimetre-wave 
counts and background intensities are required in order to 
investigate the history of star formation at earlier epochs. 



3.4 Fitting the evolution of distant galaxies 

In order to impose limits to the form of evolution at high red- 
shifts, the peak and anvil models of pure luminosity evolu- 
tion arc used to predict both the spectrum of submillimetre- 
wave background radiation and the counts of dusty galaxies 
at wavelengths of 2.8 mm, 850 and 450 fj,m. The density of 
metals generated in high-mass stars by the present epoch is 
also calculated. The models used to derive these quantities 
are discussed by Blain & Longair (1993a, 1996). 

The results of these calculations are presented in Fig. 6 
for illustration in the anvil model, assuming the values of p 
and Td derived above, as a function of Zmax and za, the re- 
maining free parameters in the evolution function (equation 
3). Likelihood contours derived for independent fits to the 
background intensities reported by Fixsen et al. (1998), and 
a combination of the results of Schlegel et al. (1998) and 
Hauscr et al. (1998) arc compared with the likelihood of a 
fit to SIB's 850-/^m counts. 

Additional constraints can be imposed by considering 
source counts at wavelengths of 2.8 mm (Wilner & Wright 
1997) and 450 /xm (SIB; Barger et al. 1998), and the den- 
sity parameter of heavy elements generated by high-mass 
stars by the present epoch f2in(0), expected to be approxi- 
mately 10~^. If solar metallicity (Savage & Sembach 1996) 
is typical of the Universe as a whole, and the density param- 



eter in baryons fib/i^ = 0.019 (Buries & Tytler 1998), then 
flm — 1.9 X 10^^. These constraints arc weaker than those 
imposed by the counts and background data and shown in 
Fig. 6. The background estimate of Puget et al. (1996) pro- 
vides a substantially weaker constraint on our models than 
any of the other estimates of the background and so does 
not add anything to our analysis. As noted by Fixsen ct al. 
(1998), if the H"*" correction employed by Puget et al. (1996) 
is omitted, then the two determinations of the background 
intensity are in good agreement (see Fig. 4b of Fixsen et al. 
1998). The results of further analysis of the COBE data will 
probably increase the certainty with which the background 
is known, but, £is shown in Fig. 6 and discussed below, the 
newly determined submillimetre-wave counts are the crucial 
data for discriminating between models. 

In Fig. 7 the joint likelihood for a fit to all the data, 
apart from Puget et al.'s background is presented for both 
the peak and anvil models. All the data can be fitted very 
satisfactorily. Note that changing the values of the param- 
eters p and Td, which were determined using low-redshift 
data in Section 3.3, within their specified errors causes the 
probability contours in 2max-2o parameter space to move by 
less than Icr from the positions shown in Fig. 7. 

There are several interesting points to note in Fig. 6. 
First, the regions of parameter space in which both the in- 
tensity of background radiation and the 850-/im counts are 
fitted overlap satisfactorily. The counts and background es- 
timates are consistent, although models that fit the counts 
arc in better agreement with the recent background esti- 
mates by Fixsen et al. (1998) than with Puget et al. (1996). 
Secondly, the tightness of the constraints imposed on the 
model parameters by the background measurements and the 
850-/im counts are comparable. As the accuracy of the 850- 
/im counts is set to increase significantly when the complete 
cluster sample is presented (Blain et al. submitted; Small et 
al. in preparation) and the data firom other surveys becomes 
available, knowledge derived from the counts will be the de- 
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Solid: anvil model 
Dashed: peak model 
Dotted: no evolution 



-1 -0.5 0.5 1 

log^Q(Fhix Density / Jy) 

Figure 5. Observed 60- //m counts of /i?j4 5 galaxies and the best- 
fitting peak and anvil models determined from the analysis in 
Section 3.3, plotted in the format used by Oliver et al. (1992). 
For compajrison, the counts expected in a non-evolving model are 
shown by the dotted line. Tiic data arc taken from Hacking & 
Houck (1987; crosses), Rowaii-Hobinson ct al. (1990; empty tri- 
angles), Saunders et al. (1990; filled triangles) and Gregorich et 
al. (1995; circles): see also Berlin et al. (1997). 




Figure 6. The relative significance of the high-redsliift con- 
straints imposed on models of galaxy evolution that account for 
the low-redshift data. The results are derived in the anvil model 
for illustration. The corresponding diagram in the peak model 
has a very similar form. The results of fitting the FIRAS back- 
ground radiation intensity determined by Fixsen et al. (1998), the 
DIRBE background (Schlegel et al. 1998; Hauser et al. 1998) and 
SIB's 850-/im counts are compared. The counts are represented 
by the extremely thick solid line. The probability of fitting the 
FIRAS background is represented by the shaded regions; a mini- 
mum reduced value of 0.07 is obtained for the best fit to this 
background. The probability of fitting the DIRBE background is 
represented by the tramlines. 



termining factor in discriminating between models of galaxy 
evolution. Thirdly, the locus of the best-fitting parameters 
is a line and not a distinct peak. Hence, given the existing 
data, no individual model provides a best fit, but rather fam- 
ilies of acceptable models trace out this locus (see Fig. 7) . In 



the final section we discuss how we can discriminate between 
members of these families. 

Three values of zo and Zmax are selected in each 
model, highlighted by dots in Fig. 7. The models defined 
by these sets of parameters arc consistent with all the cur- 
rent submillimetre-wave data, and we list their parameters 
in Table 1. 

Submillimetre-wave observations are very sensitive to 
the form of galaxy evolution at high redshifts, and so provide 
a scries of interesting cosmological tests (Blain 1998a). At 
present we are not prepared to use the limited available data 
to infer values of Qo and Qa; however, note that low values 
of Clo or a non-zero value of tend to increase the back- 
ground radiation intensity for a given count. In an Einstein 
de Sitter model, an exact fit to the counts requires a more 
significant form of evolution than an exact fit to the back- 
ground (see Fig. 6). Hence, at a first level of approximation, 
a high-density Universe would seem to be more favoured by 
this data. We will return to this issue once more accurate 
counts have been determined. 



3.5 An optically-motivated model 

The likelihood contours in Fig. 7 have been used to select 
six best-fitting peak and anvil models, to test the robust- 
ness of the procedure, and to investigate whether or not 
observational tests can be devised to distinguish between 
the different cases. However, these best-fitting models have 
been selected on the basis of observations in the far-infrared 
and submillimetre wavebands alone. If the extinction cor- 
rections applied to account for reddening in faint Lyman- 
dropout selected sources (Pettini et al. 1998a,b) are correct, 
then the empty squares in Fig. 1 represent the global SFR 
in a high-redshift population of non-active sources, with in- 
dividual SFRs in the range 10 to SOMoyr"^ The effects 
of this population, which consists of less luminous sources 
than the galaxies detected by SIB, should be included in the 
models. It is possible to use a peak model - Peak-G in Ta- 
ble 1 - to connect this population and the low-redshift IRAS 
galaxies, assuming that they are subject to the same form 
of low-redshift evolution as the IRAS galaxies. The choice 
of zq in this model is unimportant, as the SFR at > 5 is 
effectively zero in this model. 

In order to explain the high-redshift data - the 
submillimetre-wave background radiation intensity and 
counts, and the majority of the 175-/im ISO counts - an ad- 
ditional second population of more luminous, but obscured, 
galaxies must be included. This 'Gaussian' population is 
based on the IRAS luminosity function, and is assumed to 
undergo pure luminosity evolution of the form. 



gg{z) = eexp<^ - 



[t(z)_t{z^ 
20-2 



(5) 



This function represents a burst in the SFR of the pop- 
ulation at a characteristic epoch, centred on the cosmic 
epoch that corresponds to redshift Zp, with a characteris- 
tic timescale a and a normalization with respect to the 
low-redshift IRAS luminosity function. To avoid introduc- 
ing extra parameters, the dust temperature in this addi- 
tional population remains at 38 K. The three parameters 
0, Zp and CT are determined by fitting the model to the 
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submillimetre-wave counts and background radiation inten- 
sity (see Table 1). Reasonably, Qa, which is proportional to 
the total amount of star formation activity in the burst, is 
constrained more tightly than the individual values of and 
a: Qa = 76±8Gyr. 



4 CONSISTENCY ARGUMENTS 

4.1 Generation of dust and stars 

Wc have adopted a singlo-temperaturo dust model to eval- 
uate a star formation history in obscured galaxies that is 
consistent with all existing observations at wavelengths be- 
tween the millimetre and mid-infrared wavebands. The tem- 
perature was obtained in Section 3.3 using the low-redshift 
count data. For these models to be self-consistent, a suffi- 
ciently large mass of dust must be present at each epoch 
to absorb the optical/ultraviolet continuum luminosity of 
young star-forming regions and re-emit it at the required 
temperature. This type of condition was discussed in the 
context of a starburst in an individual galaxy by Eales & Ed- 
munds (1996, 1997). While preparing previous papers (Blain 
& Longair 1993a, 1996), the results of 'back-of-the-envelope' 
calculations suggested that this effect would not be ijartic- 
ularly significant, especially as there was so little data with 
which to constrain model predictions. However, now that 
much more data is available, it is important to confirm the 
self-consistency of our models. 

The relationship between the mass, temperature and 
total bolometric luminosity of dust grains is reasonably well 
determined. If the quantities are all expressed in densities 
per unit comoving volume, then the luminosity density emit- 
ted in the rest-frame far-infrared waveband 

( ^ :^ 2 X 10-^ ( ^^^"1 ( ^) , (6) 

where pd is the comoving density of dust grains, and a u^'^ 
cmissivity law is assumed. Regardless of whether the dust 
is distributed uniformly or in clumps, equation (6) implies 
that there is a well-defined minimum temperature that the 
dust must have in order to emit a given luminosity density. 
Because of the strong temperature dependence on the right- 
hand side, it is possible for a small fraction of the dust at a 
higher than average temperature to dominate the emission 
of energy; for example, the same luminosity is emitted by a 
mass of dust at a uniform temperature as by only 1 per cent 
of the grains at a teniperaturc increased by a factor of 2.3. 
However, it is not possible for a required luminosity to be 
generated at a uniform lower temperature. 

For an obscured star formation model to be self- 
consistent, enough dust must be generated by prior genera- 
tions of stars to satisfy the constraint on dust temperature 
imposed by equation (6) at a temperature of 38 K. If this 
condition is not satisfied, then insufficient absorption can 
take place at high redshifts. By rearranging equation (6), 
the limit is given by 

/ \ o / \ -0.182 

(^)>7.i( ^^'^ 3 ) ( — e^] . (7) 

Pfir is expected to be directly proportional to the SFR 
per unit comoving volume psfr; thus 



pFIR — KpsFB.- (8) 

Thronson & Telesco (1986) derived a value of 7^ = 1.5 x 
10"' L0 Mq^ yr for the entire range of stellar masses included 
in a Salpeter initial mass function (IMF). We note that the 
choice of this IMF is somewhat arbitrary at this time owing 
to our lack of knowledge of the IMF in galaxies in the early 
Universe. Rowan-Robinson et al. (1997) suggest a value of 
K = 2.2 X 10^ L0 Mg^ yr, and we assume this value here. If 
only OBA stars are formed, then the value of K is expected 
to increase by a factor of about 3. Equations (6) and (8) 
provide a simple relationship between the bolometric lumi- 
nosity emitted by dust and the comoving density of stars 
formed. Hence, given any star formation history, the total 
density of material converted into stars by any epoch 

rt(zo) 

p'{z) = K~'^l pFmdt-l-piii, (9) 

Jt{z) 

can be calculated, pjjj is the comoving density of stars 
formed in any Population HI component at redshifts greater 
than zo. 

The evolution of the density of dust, which is required 
to evaluate the limit in equation (7), can also be estimated 
using a simple model. Throughout the process of galaxy for- 
mation, dust is created as stars die and redistribute metals 
into the interstellar medium (ISM) . Dust is also expected to 
be destroyed either by being bound up in newly formed stars 
or evaporated by the intense ultraviolet/blue light from OB 
stars. In general, both processes would be expected to pro- 
ceed at a rate approximately following the SFR, and so it 
seems reasonable to assume that 

Pi(z) ~ £dP*(^), (10) 

in which is the efficiency of dust formation: that is, if a 
unit mass of gas is converted into stars, then a fraction td of 
this mass is returned to the ISM in the form of dust grains 
after a short time. The details of the processing of elements 
in stars and the retention of metals in stellar remnants are 
extremely complex, but an efficiency of between 10"'' and 
10~^ seems reasonable. If we assume that 30 per cent of the 
material in stars with masses greater than 20 M© is con- 
verted into metals and returned to the ISM, then for either 
a Salpeter IMF, with mass limits of 0.07 and 100 M©, or a 
Kennicutt IMF (Kennicutt 1983), with mass limits of 0.1 
and 100 Mq, 2 per cent of the material processed into stars 
is returned to the ISM in the form of metals. Note that 
changing the density of dust in equation (7) by an order of 
magnitude modifies the lower limit to the dust temperature 
by only a factor of 1.5. 

Hence, by combining equations (7), (8), (9) and (10) 
a lower limit to the dust temperature that is required for a 
self-consistent model can be derived. If this lower limit is less 
than the rcdshifted temperature of the CMBR Tcmbr(1+2), 
then clearly the limit is not a practical restriction (Eales & 
Edmunds 1996). If the lower limit to the dust temperature 
is calculated assuming that the dust grains are in thermal 
equilibrium with the CMBR, then a dust temperature of 
38 K at low redshifts would correspond to a slightly higher 
temperature of 40 K at z = 10. 
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Figure 7. The joint probability that parameters defining the anvil and peak models fit to all the high-redshift data listed in 
Section 3.4, apart from the background intensity of Puget et al. (1996). The maximum-likelihood values, marked by stars, correspond 
to values of 1.15 and 1.12 in (a) and (b) respectively. The chosen parameters in the models used in subsequent sections (Table 1) 
are indicated by the labelled dots. . 



Table 1 . The parameters defining a series of models that adequately describe the star formation histories derived from optical surveys and 
from the submillimetre-derived likelihood contours shown in Fig. 7. The Peak-M and Peak-P models fit the optically-derived histories: 
see Fig. 1. No extinction correction is included in the Peak-M model (Madau et al. 1996). The Peak-P model includes an extinction 

correction (Pcttini ct al. 1998). The Peak-G model is chosen to fit both the low-rcdsliift evolution and normalisation of the luminosity 
density of IRAS galaxies and the extinction-corrected optical points (Pettini et al. 1998a, b): see Section 3.5. The other models fit all the 
submillimetre/far-infrared data. The line styles and thicknesses listed refer to the form in which the results of these models are presented 
in Figs 8 to 13. = 38 ± 4 K in all models. 



Model p Hm / tsfr(O) / 10 ^ Line style Line Comments 

10~3 MQyr-iMpc~3 Thickness 



Peak-M 


10 


0.39 ±0.02 


4.6 ±0.3 


0.16 


8.0 ±0.8 


Dot-dashed 


Thin 


Peak-P 


10 


0.48 ± 0.05 


4.9 ±0.4 


0.31 


7.2 ±0.8 


Dot-dashed 


Thick 


Anvil-10 


10 


1.49 ±0.15 


3.8 ±0.2 


0.85 


7.7 


Solid 


Thick 


Anvil-5 


5 


1.81 ±0.15 


3.8 ±0.2 


0.99 


7.7 


Solid 


Medium 


Anvil-4 


4 


2.03 ±0.15 


3.8 ±0.2 


1.05 


7.7 


SoHd 


Thin 


Peak- 10 


10 


1.04 ±0.15 


3.9 ±0.2 


0.79 


7.7 


Dashed 


Thick 


Peak-5 


5 


1.17±0.15 


3.9 ±0.2 


0.84 


7.7 


Dashed 


Medium 


Peak-4 


4 


1.42±0.15 


3.9 ± 0.2 


0.90 


7.7 


Dashed 


Thin 


Peak-G 


10 


0.63 ±0.04 


3.9 ±0.1 


0.35 


7.7 


Dotted 


Thin 


Gaussian 


10 


N/A 


N/A 


0.91 


N/A 


Dotted 


Thick 



Fits Madau et al. (1996) 
Fits Pettini et al. (1998a,b) 



Base for Gaussian model 

Always added to above 
Zp = 2.1 ± 0.2 
a = 0.8±0.5Gyr 
e = 95 ± 40 



4.2 The temperature limit for selected models 

The redshift-dependence of the dust temperature limit cal- 
culated for all seven models listed in Table 1 that axe fitted 
to the submillimetre-wave data is shown in Fig. 8. In each 
case Ed ~ 10~^. In Fig. 8(a) pin = 0. With no Population III 
stars to generate early dust, the lower limit to the tempera- 
ture exceeds the assumed value of 38 K at redshifts close to 
zo in all the models. The consistency breats down only at 
redshifts greater than about 8 for the models with zo = 10. 
Only a small fraction - less than 10 per cent of sources 
in these models are at redshifts greater than 8, as shown 
by the redshift distributions in Fig. 13. Hence, even this vi- 
olation of consistency does not invalidate the conclusions of 



the paper. In Fig. 8(b) a small Population III is included, 
in which pjjj = 10~^p*(0). By seeding the high-redshift gas 
with heavy elements and dust, this population ensures that 
all the models are self-consistent. 

The self-consistency condition is significant only if a low 
typical temperature of dust, such as 38 K, is assumed. If a 
significant fraction of submillimetre-wave sources are heated 
by AGN rather than high-mass stars, then systematically 
less dust forms by any epoch, and so the condition would be 
more difficult to satisfy. In an hierarchical model, a certain 
fraction of star formation is likely to occur in brief bursts 
at the epochs of mergers, and so only a fraction of the total 
mass of dust will be emitting strongly at any epoch. The 
lower limit to the dust temperature imposed by equation 
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Temperature in model 














Figure 8. The lower limits to the dust temperature imposed as a function of redshift by the consistency argument detailed in 
Section 4. The temperature of 38 K assumed in the models and the lower limit imposed by the CMBR temperature Tcmbr a-re 
also shown. In a self-consistent model, the lower limit to the temperature must be less than 38 K at all epochs. The line styles 
and thicknesses axe interpreted in Table 1. In (a) no Population-Ill stars are included. In (b) an additional 2 per cent of the total 
amount of star formation taJiing place in each model is added to generate dust in a Population III at redshifts greater than 10. 
When this additional small population is included, the models are all self-consistent. 



(7) would also increase in this case. In the models of Baugh 
et al. (1998) less than 20 per cent of stars form in such bursts, 
but the bursts could be responsible for a larger percentage 
of submillimetre-wave sources. Either a large population of 
dusty AGN or hierarchical mergers at high redshifts would 
make it more difficult to obtain a self-consistent model. 

Based on the calculations discussed above all the mod- 
els listed in Table 1 are self-consistent if a few per cent of 
star-formation activity occurs in Population III stars. If this 
population does not exist, then the first generation of stars 
would not be obscured by dust for a brief period. Neverthe- 
less, even in this case the results in this paper would not be 
modified by more than about 10 per cent. 



5 RESULTS 

5.1 Introduction 

Here we use the models presented in Section 3 that account 
for all the available data in the millimetre/submillimetre 
and far-infrared wavebands to estimate the star formation 
history in obscured sources, and relate it to the observed 
consumption rate of neutral gas (Storrie-Lombardi, McMa- 
hon & Irwin 1996) and metal enrichment (Pettini et al. 1997) 
in Lyman-a absorbers. 

5.2 The star formation history 

The star formation histories associated with the models pre- 
sented in Table 1 are plotted in Fig. 9. Although there is a 
reasonable scatter in the results that fit the submillimetre- 
wave data, all the predictions lie within a factor of about 3 
of each other over the critical range of redshifts from 2 to 
5. The offset between the optically derived points and the 
submillimetre-derived curves is uncertain to within a fac- 
tor of about 2. The estimated uncertainties in the values of 
the fitted parameters listed in Table 1 can shift the curves 
through a region similar to that enveloped by the curves 



from model to model. The histories are fully consistent with 
limits derived from existing background radiation intensities 
(Blain & Longair 1993b; Burigana et al. 1997). 

At low redshifts the average slope of the stax formation 
histories derived here is in good agreement with the form 
determined from observations of optically selected galaxies. 
However, there is a clear discrepancy between the two classes 
of history at high redshifts. The histories cannot be recon- 
ciled, even if the properties of the high-redshift, optically 
selected galaxies are corrected for the estimated dust extinc- 
tion (Pettini et al. 1998a,b; Calzetti et al. 1995). A further 
correction to the dust extinction in the ultraviolet-selected 
samples appears unlikely, and we therefore suggest that the 
remaining star formation arises in a population of obscured 
galaxies which are either too red or too faint in the rest 
frame ultraviolet waveband to be selected by the Lyman- 
dropout technique (Steidel et al. 1996b). This suggestion 
can be tested by carrying out optical and near-infrared spec- 
troscopy of submillimetre-selected sources to determine the 
fraction that lies at high redshifts 2^3. 

The limited experience so far gained from follow-up 
observations of sources selected on the grounds of their 
powerful dust emission at more modest redshifts, IRAS 
F 10214+4724 (Rowan-Robinson et al. 1991; Close et al. 
1995) and SMM J02399-0136 (SIB; Ivison et al. 1998b), 
suggests that such sources can be detected and studied 
in the optical and near-infrared waveband in reasonable 
integration times. Hence the prospects for following up 
submillimetre-wave surveys and determining the typical 
properties of the optical counterparts appear to be quite 
good. Direct submillimetre-wave observations of the Lyman- 
break objects on which the optically derived star formation 
history is based will soon allow us to determine whether or 
not the extinction corrections inferred from their optical and 
near- infrared spectra are correct, and thus the contribution 
of such sources to the population of submillimetre-selected 
objects. 

Another population which may provide a link between 
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Figure 9. Star formation histories in optically-selected and obscured galaxies presented for the models listed in Table 1. The data 
points are identical to those plotted in Fig. 3. The seven upper curves represent star formation histories that are consistent with all 
the available far-infrared and submillimetre-wave data; the three lower curves represent star formation histories that are consistent 
with the optical/ultraviolet /near-infrared data. 



the optical and submillimetre views of the distant Universe 
are the class of extremely red objects (EROs). These have 
been detected in deep multi- waveband surveys (Graham & 
Dey 1996; Ivison et al. 1997; Cimatti at al. 1998; Dey et 
al. 1998), but not in sufficient numbers to be included in 
the derivation of the points in Fig. 9. Nevertheless, detailed 
observations of these strongly obscured galaxies may help to 
unify the two regimes. 

5.3 Enrichment 

Observations of Lyman-Q absorbers along the line of sight 
to distant quasars allow the evolution of the mass of neutral 
gas and the typical metallicity in the Universe to be traced 
as a function of epoch (Storrie-Lombardi et al. 1996; Pettini 
et al. 1997). 

5.3.1 Gas consumption and the density parameter in stars 

The evolution of the neutral gas content of the Universe has 
been determined by Storrie-Lombardi ct al. (1996). The co- 
moving density parameter of neutral gas fig appears to peak 
at a redshift of about 3. When combined with an estimate of 
the density parameter in the form of stars 0,(z), measured 
at the present epoch as (0) = (5.9 ± 2.3) x 10"^ (Gnedin 
& Ostriker 1992), this result can be viewed in the context of 



the star formation history. In Fig. 10(a) the mass of material 
that has been processed into stars is derived as a function 
of epoch using equation (9) for each of the star formation 
history models listed in Table 1. 

The optically derived star formation history, as pre- 
sented by Madau ct al. (1996, 1998) appears to be very 
consistent with the data presented in Fig. 10(a). With or 
without dust corrections, that is in the Peak-M or Peak-P 
model respectively, $7, (z) builds up smoothly and O* -|- Og 
remains approximately constant as a function of epoch to 
within the errors. The same is true for the Peak-G model. 
However, the values of f2, (0) predicted by the submillimetre- 
derived star formation histories are considerably larger. The 
total density parameter of processed material exceeds the 
stellar density derived by Gnedin & Ostriker (1992) by a 
factor of about 5. In fact, the values of O*(0) determined 
in the submillimetre-derived models are comparable to the 
value of Qh. Note, however, that despite this large predicted 
stellar mass, about half of the total mass is expected to 
have been processed prior to a redshift of 2 to 3 in all the 
submillimctrc-dcrivcd histories. The observed turn-over in 
the neutral gas fraction and the maximum rate of star for- 
mation in these histories are coincident in redshift. 

If the background radiation intensity and counts deter- 
mined in the submillimetre waveband are not in error, then 
either a large fraction of the baryons in the Universe have 
been processed into stars, or the process of star formation 
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in the distant sample must be different from that in nearby 
galaxies, with an initial IMF biased to high-mass stars. A 
third possibility is that some fraction of the luminosity of 
distant dusty galaxies could be produced by AGN rather 
than by high-mass stars. 

First, let us consider the effects of modifying the IMF. 
If a Salpeter IMF is assumed, with a lower mass limit of 
0.07 Mq, then in all the models listed in Table 1 about 65- 
70 per cent of all stars formed are still burning at the present 
epoch. In the submillimetre-derivcd models this fraction is 
typically 2-3 per cent smaller than in the optically-derived 
models, certainly not by enough to account for the factor of 
about 4 by which Q, exceeds the observed value (Gnedin & 
Ostriker 1992) . If a shallower low-mass IMF is assumed, for 
example Kennicutt (1983), the fraction of stars still burn- 
ing is reduced to about 40 per cent. These results are un- 
affected by the value of the upper mass limit in the IMF. 
A population of cxtragalactic stars, not detected by Gnedin 
& Ostriker, could account for some of the difference. In or- 
der to account for a factor of 2.5 discrepancy between the 
submillimetre-derived star formation histories and Gnedin & 
Ostriker (1992) that remains after the fraction of 30 per cent 
of stellar remnants is removed from the curves in Fig. 10(a), 
a lower mass limit of about 0.7 M© is required in the IMF. 
This is certainly a plausible scenario: a modified IMF in dis- 
tant star-forming galaxies could account for these results, 
and would predict a stellar density at the prcscut epoch for 
the submillimetre-derived models that is smaller by a fac- 
tor of 4 compared with that shown in Fig. 10(a). A lower 
mass limit to the IMF of 3 M0 was suggested by Zepf & Silk 
(1996) in order to explain mass-to-light ratios of elliptical 
galaxies, and Rieke et al. (1993) interpret observations of 
M82 using the same cutoff value. The potential variations 
in the high-redshift IMF are discussed by Larson (1998). 

Alternatively, we consider the effects of a population of 
obscured sources powered by AGN rather than high-mass 
stars. Such a population almost certainly exists. Both the 
brightest submillimctre-wave object detected by SIB and 
the prototype of ultraluminous distant dusty galaxies, IRAS 
F 10214-1-4724, contain AGN which must contribute to the 
heating of dust in these galaxies (Ivison et al. 1998b). In 
order to match the observed stellar density at the present 
epoch to that predicted in the models, but without modify- 
ing a standard Salpeter IMF, about 60 per cent of the energy 
emitted by distant submillimetre-selected sources must be 
generated by AGN. Such a large fraction would require the 
lower limits to the dust temperatures derived in Section 4 to 
be increased by a factor of about 20 per cent, but this is not 
unreasonable. Sanders & Mirabel (1996) report that about 
a third of ultraluminous galaxies detected by IRAS appear 
to contain an AGN. 

The relative importance of these alternative explana- 
tions cannot be determined at present. Optical follow-up 
observations of submillimetre-selected sources will help to 
discriminate between the high-mass IMF and AGN mod- 
els. The application of more advanced chemical enrichment 
models will also provide useful information when follow-up 
observations are available (Pei & Fall 1995; Bales & Ed- 
munds 1996, 1997). 



5.3.2 Generation of metals 

The mass of metals generated in any star-formation history 
can be calculated using equations (9) and (10). If divided 
by the solar abundance, assuming that the sun formed out 
of representative gas at « = 0.38, the abundances of heavy 
elements are predicted to evolve as shown in Fig. 10(b) for 
the models listed in Table 1. The data on metal abundances 
in damped Lyman-a systems determined by Pettini et al. 
(1997) are also plotted. As they noted, the fits to the opti- 
cally derived counts and the Peak-G model are in full agree- 
ment with the data. By contrast, metals appear to be over- 
produced by a factor of about 5 at redshifts of 2 and 3 in the 
submillimetre-derived histories. However, in these histories 
the rate of enrichment shown by the curves in Fig. 10(b) is 
broadly consistent with the slope interpolated between the 
three highest redshift data points. 

Is this difference in the rate of production of metals 
significant? Clearly, modifying the low-mass cut-off of the 
IMF as discussed above will not affect the metal enrichment 
rate. However, including a more substantial population of 
AGN-powered obscured sources at high redshifts would de- 
lay metal enrichment as compared with a model in which the 
sources were powered by star formation alone. If AGN out- 
numbered starbursts by a ratio of 4:1, then the enrichment 
curves in the submillimetre-selected models in Fig. 10(b) 
would fall into line with the data. Note that this fraction 
of AGN is greater than that required to explain the discrep- 
ancy between the observed and predicted stellar density in 
Section 5.3.1. 

Alternative arguments could also be used to explain 
at least some of the excess metal density predicted in the 
submillimetre-selected models. The observations of damped 
Lyman-a systems, from which the points in Fig. 10(b) are 
derived, may sample the outer regions of disk systems (Fer- 
guson, Gallagher & Wyse 1998), rather than the highly- 
enriched cores, and so the metallicity of these systems may 
be systematically undersampled. In addition, quasars whose 
line of sight is intersected by the most metal-enriched sys- 
tems arc more likely to be obscured and thus less likely to be 
included in optically selected samples (Pei & Fall 1995). If 
this were true, then the data in Fig. 10(b) would tend to be 
obtained in directions with a systematically low metal den- 
sity, again under-representing the true smoothed fraction of 
metals. The metal abundance in intracluster gas observed by 
the ASCA X-ray observatory is typically about a third of the 
solar abundance (Gibson, Loewenstein & Mushotzky 1998), 
a value which lies above the points in Fig. 10(b), perhaps 
indicating that the global metal abundance at redshifts of 2 
to 3 lies closer to the values predicted by the submillimetre- 
derived models. 



5.3.3 Enrichment discussion 

The star formation histories in dusty galaxies that account 
for all of the far-infrared and submillirnctrc-wavc data ap- 
pear to be inconsistent with observations of (0) and high- 
redshift metal enrichment. The results depend on the as- 
sumed IMF in dusty galaxies, and the optically derived data 
with which they are compared could well be subject to sys- 
tematic biases. These biases would tend to lower the mear 
sured enrichment at high redshifts, and lead to a global value 
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Figure 10. (a) The density parameter of gas and stars as a function of redshift. The curves represent the density parameter of 
material processed into stars in the models listed in Tabid, assuming a Salpeter IMF with mass limits of 0.07 and 100 Mq. The 
data for the total density of stars at the present epoch (shaded region) were obtained by Gnedin Sz Ostriker (1992). The data 
points represent the density parameter in neutral hydrogen (Storrie-Lombardi et al. 1996). The lower three curves correspond to 
the models that fit the optically derived star formation histories listed in Table 1. (b) The rate of increase of metallicity expected 
in the same models. The data points were determined by Pettini et al. (1997). 



of enrichment in better agreement with the predictions of 
the submillimetre-derived models, which we expect to be 

unaffected by cither systematic biases or local effects. 

If the optical data is taken at face value, then either 
an IMF biased to high-mass star formation or a significant 
fraction, perhaps up to 60 per cent, of distant ultraluminous 
submillirnetre-wave sources must be powered by AGN in or- 
der to reconcile the submillimetre-wave and optical data. If 
the fraction of AGN at high redshifts reached about 60 per 
cent, then the extinction-corrected star-formation history 
discussed by Pettini et al. (1998a, b) would account for most 
of the star-formation activity that has taken place in the 
Universe. 

More detailed follow-up observations of submillimetre- 
selected sources and detailed stellar and chemical evolution 
studies will be required in order to disentangle these effects. 
In particular the fraction of AGN uncovered in follow-up 
observations of submillimetre-wave surveys will provide a 
powerful test of these ideas (Ivison et al. 1998b). 



5.4 Contamination by AGN 

Several circumstantial arguments can be used to illuminate 
the fraction of the SCUBA galaxies in which the dust is 
heated by AGN as compared with star formation using ex- 
isting data. 

The proportion of AGN that are missing from optically 
selected catalogues due to obscuration by dust has been dis- 
cussed by Fall & Pei (1993). The absorption of the optical 
emission from AGN by intervening dust is also discussed by 
Dririkwater et al. (1996). Fall & Pel's models suggest that 
the true surface density of quasars with luminosities greater 
than of order 10^^ L© at redshifts in the range 2 to 3 - suf- 
ficiently luminous to be detected by SCUBA - should be of 
order 50 deg^^, after correcting for the effects of dust extinc- 
tion. For comparison the surface density of SCUBA galaxies 
is about 2000 deg"^. 

Based on a radio selection technique, insensitive to the 



effects of dust absorption, Webster et al. (1995) estimate 
that 80 per cent of radio-loud quasars are missing from 
optically selected samples. If the same fraction of radio- 
quiet quEisars is absent, then the surface density of obscured 
quasars that are sufficiently bright to be detected by SCUBA 
should be of order lOOdeg"^. The results of deep X-ray sur- 
veys are also insensitive to the effects of obscuration by dust, 
and resolve a significant fraction of the X-ray background ra- 
diation. In order to account for the background spectrum, a 
population of dust obscured AGN, which are not detected in 
optical surveys, must be included in models (Almaini et al. 
1997, 1998; Almaini, Lawrence & Boyle 1999). The surface 
density of faint X-ray sources, which are candidates for this 
population, is about 230±40deg~^ (Jones et al. 1997). 

These observations all tend to suggest that the surface 
density of luminous dust obscured AGN might be typically 
up to a few tenths of the surface density of SCUBA galaxies, 
and thus that star-formation activity is responsible for the 
greater fraction of the submillimetre-wave flux density emit- 
ted by the population of SCUBA galaxies. This proportion 
is consistent with the fraction of compact sources identified 
with SCUBA galaxies by SmaU et al. (1998). 

Only one SCUBA-selected source, SMM J02399-0136, 
has been subject to detailed follow-up observations (Ivison 
et al. 1998). While SMM J02399-0136 certainly contains an 
AGN, a significant fraction of its luminosity is thought to be 
produced by massive star formation. The recent detection 
of redshifted CO(3^2) emission from SMM J02399-0136 
using the Owens Valley Millimeter Array (Frayer et al. 1998) 
supports the idea that approximately equal amounts of its 
luminosity are contributed by star formation activity and an 
AGN. 

If the SCUBA sources are powered predominantly by 
star formation activity, then a global gas consumption rate 
of order 0.5 M© yr^"*^ Mpc^'^ is required at redshifts of about 
3. This rate can be compared with the predictions of the 
self-consistent chemical evolution models presented by Pei 
& FaU (1995). The infall model of Pei & FaU (1995) predicts 
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Figure 11. Predicted background radiation intensities in the 
eight models of galaxy formation listed in Table 1 that are fit- 
ted to submillimetre-wave data. The background data is plotted 
in the format of Fig. 2. 

a very similar consumption rate, and so there is no practical 
limit to the fraction of the luminosity of SCUBA sources 
that can be generated by star formation rather than AGN. 
Detailed multiwavebarid follow-up observations will be re- 
quired to determine the relative importance of AGN and star 
formation for heating dust in SCUBA galaxies. At present 
it seems likely that up to about a third of the emission from 
these sources is powered by AGN, with the remainder due 
to star formation. Interestingly, this is the fraction of the 
count of SCUBA galaxies predicted on the basis of the X- 
ray background and luminosity function by Almaini et al. 
(1999). 



6 THE CONSEQUENCES FOR FUTURE 
OBSERVATIONS 

The submillimetre-derived star formation histories shown in 
Fig. 9 correspond to models of galaxy evolution that ac- 
count for the low-redshift evolution of IRAS galaxies, the 
observed submillimetre-wave background radiation inten- 
sity and the 850-/im counts, and predict counts which are 
in agreement with the limits determined at wavelengths of 
2.8 mm and 450 nm. The predicted background radiation in- 
tensities, source counts and redshift distributions at a wide 
range of wavelengths are derived in these models and the 
results are used to discuss the future of submillimetre-wave 
cosmology. The most promising observations with which to 
obtain more accurate determinations of the history of ob- 
scured star formation are considered. 

6.1 Background radiation 

Predictions for the intensity of background radiation in the 
millimetre/submillimetre and far-infrared wavebands de- 
rived in the models listed in Table 1 are shown in Fig. 11. 
Similar spectra are predicted by all the models, regardless 
of their very different star formation histories. The ensemble 



Table 2. The redshifts below which 10, 50 and 90 per cent of 
sources detected in surveys at reasonable flux density limits at a 
range of wavelengths are expected to lie, zio , 250 and zgo respec- 
tively. The values are obtained from the spread of the redshift 

dis! ribnt ioi) curA'c^s in Fi,;^;. 13. 
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of models included in Fig. 11 predicts background intensities 
of (8.3 ± 1.2) X 10"^° and (4.6 ± 0.8) x 10"^ Wm'^ sr'^ at 
wavelengths of 850 and 450 ^m respectively. For comparison, 
the corresponding values reported by Puget et al. (1996) and 
Fixsen et al. (1998) are (2.7±2.0) x 10"^°, (2.3±2.0) x 10"^ 
and (5.3±1.6) x 10"^°, (3.3±1.0) x 10"^ Wm'^ sr'^ respec- 
tively. The form of the background spectrum is therefore not 
a good quantity with which to discriminate between differ- 
ent models of the star-formation history. 

6.2 Source counts 

The source counts predicted at wavelengths of 2.8 mm, 

1.3 mm, 850 /im, 450 /im, 175 /im, 60 /im and 15 /im are 
shown in Fig. 12. The observed counts derived by Wilner 
& Wright (1997), SIB, Holland et al. (1998b), Kawara et 
al. (1997), the authors listed in Fig. 5, and Oliver et al. 
(1997) arc also plotted. The predicted abundance of lensed 
submillimetre-wave galaxies (Blain 1996) is discussed else- 
where (Blain 1998c). 

The 15-//m counts are included for completeness only. 
The models used here have little predictive power at wave- 
lengths shorter than about 40 ^m, and correspondingly the 
15-^m counts have little predictive power for determining 
the star formation history at high redshifts. The median 
redshift of the sources detected by ISO at 15 /xm is expected 
to be less than 0.5 (see Fig. 13a). The 15-/im counts are uni- 
formly underpredicted by a factor of about 3, about the size 
of the observation uncertainty. The predicted count curves 
can be made to traverse the 15-/um count error bar shown 
in Fig. 12(b), by sweeping the value of the high-frequency 
power-law index in the spectral energy distribution model 
(see Section 3.2) through the range —1.4 to —2.2. The value 
of this power-law index has almost no effect on either the 
source counts or background radiation spectra expected at 
longer wavelengths. 

The spread in the counts predicted across the range 
of well-fitting models listed in Table 1 offers excellent 
prospects for discriminating between the models using 
submillimetre-wave counts determined at a range of flux 
densities and observing wavelengths. A range of exist- 
ing and future facilities will be useful for determining 
counts; Planck Surveyor (Bersanelli et al. 1996), large 
ground-based millimetre/submillimetre- wave interferometer 
arrays (Brown et al. 1996; Downes et al. 1996), advanced- 
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Figure 12. Counts predicted by the models listed in Table 1, as compared with available data. The 60-/im counts are shown in 
(a), those at 15 and 175 /im are shown in (b), those at 850 /im and 2.8 mm are shown in (c), and those at 450 /im and 1.3 mm are 
shown in (d). The references to the data in (a) are given in the caption of Fig. 5; in (b), (c) and (d) they are written adjacent to 
the data points. 



technology single-antenna ground-based telescopes (Schlo- 
erb 1997; Stark et al. 1998) and FIRST (Pibratt 1997). 

6.3 Redshift distributions 

The redshift distributions of submillimetre-selected sources 
at, or just below, the flux density limits of current surveys 
are shown in Fig. 13. The sensitivity thresholds and the red- 
shift ranges within which the redshift distributions of source 
in these surveys are expected to be 10, 50 and 90 per cent 
complete arc listed in Table 2. Based on the values listed in 
Table 2, the median redshift of SIB's submillimctrc-sclccted 
sample should be in the range 2.4 to 4.4. This is in broad 
agreement with constraints obtained from broadband pho- 
tometry of the optical counterparts by Small et al. (1998), 
as well as the first spectroscopic identifications (Ivison et al. 
1998b). 

The crucial result revealed in Table 2 is that most of the 
sources detected in submillimetre-wave surveys are likely to 
lie at redshifts well below 10, and so a large proportion will 
be accessible to multi-waveband study using 8-m-class tele- 
scopes. The advantage of a cluster lens survey, of course, 
is that in addition to the rich archival data available, the 
sources are typically magnified by a factor of a few, both 



in the submillimetre and the optical waveband. Not only 
does this make it easier to detect these faint sources in the 
submillimetre waveband, but it also facilitates spectroscopic 
follow-up observations as compared with a blank-field sur- 
vey. Blank-field submillimctrc-wavc surveys have wisely ex- 
ploited the most carefully surveyed areas of the field sky - 
the Hubble Deep Field, CFRS fields and Hawaii Deep Sur- 
vey fields (Hughes et al. 1998; Bales et al. 1998; Barger et 
al. 1998). 

If confirmed by spectroscopic observations, the rela- 
tively modest median redshift predicted for submillimetre- 
selected galaxies would further suggest that they comprise 
a distinct population of strongly obscured sources that have 
evaded detection in optical surveys on the grounds of their 
extreme obscuration rather than their extreme distance. 

6.4 Consequences for earlier results 

In earlier papers we discussed the first results of our survey 
(SIB) and the consequences of our detection of a population 
of distant dusty galaxies for source confusion in existing and 
future observations (BIS; Blain et al. 1998b). The first re- 
sults of follow-up observations in other wavebands have also 
been published (Ivison et al. 1998b). The results presented 



© 0000 RAS, MNRAS 000, 000-000 



16 A. W. Blain, Ian Smail, R. J. Ivison & J. -P. Kneib 



(a) 



Look-back time ^ Gyr 
2 4 6 8 10 12 



15 and 60/i.m ^^^^Z^" ' 






/////•' 




//// : 




''// ■•' 




iff-'- 




w 




- li 




£ 






1 0.4 and lOmJy 
r 



- 0.2 



(b) 



10 




- 0.4 



- 0.2 



(c) 



Look— back time / Gyr 
2 4 6 8 10 12 



(d) 




Look— back time / Gyr 
2 4 6 8 10 12 



6oOMm 































2mJy 



- 0.6 



0.4 



(e) 



Look-back time / Gyr 
2 4 6 8 10 12 



(f) 



1.3mm 


/ -7 ^ / 




iff / / 










■ ' ' ' / 




: '• 1 1 / 




i fi// . 




; / / / / / 
•' //' '/ / 






'j^^^t 2mJy 



- 0.8 



Look-back time / G3"r 
2 4 6 B 10 12 



2.6mm \ , 




jp— 


' ^ — f 
/ 1 

1 h 


■' / ' /' 

•' r, h 


1 I 
1 I 
1 1 


! /:/! ' 1 - 

•■ / '/' ' / 
I'll / 








2mJy 




- 0.8 



Figure 13. Predicted cumulative redshift distributions for submillimetre-selected sources detected at the wavelengths and flux 
density limits listed in Table 2. The styles and widths of the plotted curves correspond to the key in Table 1. The PeaJc-G model 
is not plotted. The redshift distributions at 15/im/0.4mjy and 60;im/10mJy are expected to be very similar. A range of surveys 
with different flux density limits at wavelengths of 450 and 850 /^m will soon be completed. The redshift distributions predicted at 
a flux density limit of 10 mjy were also calculated at both these wavelengths, and found to be very similar to those shown in (c) 
and (d). The results of such surveys can therefore be compared directly with the predictions made here. Based on the colours of 
the optical counterparts to the faint SCUBA-selected sources from the extended cluster survey described by Smail et al. (1998), 
greater than 75 ± 20 per cent and 50 ± 15 per cent of the counterparts have redshifts less than 5.5 and 4.5 respectively. 
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here reinforce the conclusions presented in these earlier pa- 
pers, demonstrating that simple well-constrained and self- 
consistent families of models of the evolution of distant ob- 
scured galaxies and quasars can account for all the available 
data in the millimetre, submillimetre and far-infrared wave- 
bands. A more complete description of the consequences for 
future surveys lies outside the scope of this paper: see Blain 
et al. (1998b) and Blain (1998c). 



7 CONCLUSIONS 

We have presented a family of self-consistent models of the 
evolution of dusty galaxies that account for all the available 
data in the submillimetre and far-infrared wavebands. 

(i) The results suggest that the luminosity density emit- 
ted in the far-infrared waveband at z > 1 evolves in a very 
different way to that inferred from observations in the op- 
tical waveband. By integrating our new luminosity densi- 
ties over all epochs, we find that the total optically-derived 
volume emissivity (Madau et al. 1996) underpredicts the 
true value by a factor of about 4. Furthermore, we can- 
not reproduce the results of submillimetre-wave observa- 
tions using a volume emissivity that evolves in the same 
way as the optically derived form. The peak of the comov- 
ing submillimetre-derived luminosity density also occurs at 
a considerably higher redshift than that suggested by optical 
observations, z ~ 2-5. 

(ii) If our submillimetre-derived star-formation histories 
are associated with a Salpeter IMF with a lower mass limit 
of 0.07 Mq , then the number of low-mass stars at the present 
epoch and the apparent metallicity at z; ~ 2-3 are overpro- 
duced by factors of up to about 5. 

(iii) In order to account for the discrepancy in the abun- 
dance of low- mass stars, then either the IMF in distant dusty 
galaxies must have a low-mass cut-off of about 0.7 M0 , or up 
to 60 per cent of the energy emitted by distant dusty galaxies 
must be generated by gravitational accretion in AGN. 

(iv) In order to understand the discrepancy between the 
metal enrichment at moderate redshifts, either gravitational 
accretion in AGN must be responsible for a large fraction of 
the energy emitted by distant dusty galaxies or the current 
optical observations are subject to considerable selection ef- 
fects that lead to the underprediction of the metal content 
at 2 ~ 2-3. 

(v) The current evidence suggests that the population of 
distant luminous submillimetre-selected sources is missing 
from surveys made in the optical waveband because of their 
very considerable dust extinction, rather than because they 
lie at extremely great distances. 

Using the results presented in this paper, we can de- 
termine the most appropriate strategy for making further 
progress in understanding the nature of distant dusty galax- 
ies and their relationship with optically-selected distant 
galaxies: 

(i) In order to learn more about the evolution of dusty 
galaxies, measurements of the submillimetre-wave counts at 
wavelengths longer than 200 ^m will allow the parameters 
that define the families of well-fitting models to be deter- 
mined more accurately. However, it will still not be easy to 



discriminate between different models within the same fam- 
ily using more accurate counts. 

(ii) In order to address the relationship between sub- 
millimetre-selected samples and optically selected sam- 
ples, follow-up optical and near- infrared imaging and spec- 
troscopy of the submillimetre-selected sources will be re- 
quired. On the basis of our predicted redshift distributions 
for the submillimetre-selected sources, it seems likely that 
a significant fraction of the sources will be detectable using 
the best existing optical and near-infrared instruments. 

(iii) The redshift distribution of submillimetre-selected 
sources will allow us to discriminate between different mem- 
bers of the same well-fitting family of galaxy formation mod- 
els that can account for all the submillimetre-wave source 
count and background radiation data. The determination of 
a substantially complete redshift distribution is the key goal 
of follow-up observations during the course of the next year. 

(iv) The optical and near-infrared spectra obtained to 
derive such a redshift distribution will be crucial for de- 
termining the relative importance of star-forming galaxies 
and AGN in submillimetre-selected samples, and therefore 
to investigate the apparent discrepancy between the rate of 
metal enrichment inferred from optical and submillimetre- 
wave observations. 

We now have a consistent picture of the evolution of the 
submillimetre-wave Universe; however, its relationship to 
our understanding derived from optical observations remains 
unclear. Some part of the discrepancy will be explained by 
the accurate quantification of known selection biases in the 
optical observations. The crucial information required to in- 
vestigate the history of star formation in more detail will be 
provided by larger submillimetre-wave surveys, and, most 
importantly, by optical and near-infrared follow-up observa- 
tions of submillimetre-selected sources. If these observations 
confirm that most submillimetre-selected sources are distant 
and show no signs of AGN activity, then there is a strong 
contradiction to investigate. 
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